Cooperative effects in the light and dark periods of two dipole-interacting atoms 
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If an atom is able to exhibit macroscopic dark periods, or electron shelving, then a driven system 
of two atoms has three types of fluorescence periods (dark, single and double intensity). We propose 
to use the average durations of these fluorescence types as a simple and easily accessible indicator 
of cooperative effects. As an example we study two dipole-interacting V systems by simulation 
techniques. We show that the durations of the two types of light periods exhibit marked separation- 
dependent oscillations and that they vary in phase with the real part of the dipole-dipole coupling 
constant. 

PACS numbers: 42.50.Ar, 42.50.Fx 



I. INTRODUCTION 



Cooperative effects in the radiative behavior of atoms stored in a trap may arise from their mutual dipole-dipole 
interaction if the atoms are close enough to each other. This is particularly interesting for two or three atoms and 
has attracted considerable interest in the literature pj- p5fl . Recently, the present authors investigated in detail the 
transition from antibunching to bunching for two two-level systems with decreasing atomic distance [|6| . 

For a single multi-level system with a metastable state one has the striking phenomenon of macroscopic dark 
periods, or electron shelving, in which the electron is essentially shelved for seconds or even minutes in a metastable 
state without photon emissions J27J- (35). For two such systems their fluorescence behavior would, without cooperative 
OS 1 effects, be just the sum of the separate photon emissions, with dark periods of both atoms, light periods of a single 
and of two atoms. In Ref. (36) the fluorescence intensity of two and three such ions in a Paul trap was measured and 
1 a large fraction of almost simultaneous jumps by two and even all three ions was recorded. This fraction was orders 
O j' of magnitudes larger than that expected for independent ions. A quantitative explanation of such a large cooperative 
! 1 „ . effect has been found to be difficult p^ , |37| -^0[| , and we will briefly discuss this question again in the last section. Other 
experiments at larger distances and with different ions showed no cooperative effects Jjj"|Jf2| . A numerical approach 
to the study of double jumps faces the problem that for good statistics one needs very long simulation times. 

As a simpler indicator of cooperative effects for systems with light and dark periods we propose here to use the 
mean durations, To, T\, and T2, of the three different types of fluorescence periods, i.e. dark, single-intensity, and 
double-intensity periods, respectively. The running time can be much shorter than required for double jumps, making 
these quantities easily accessible, both experimentally and in simulations. 

In this paper we therefore present a study of cooperative effects on the mean durations of the three types of 
fluorescence periods of two three-level V systems with a metastable state, as a function of their distance r. The level 
scheme of a single V system is depicted in Fig. 1 Q. Our simulations show that the mean durations of the single- 
and double-intensity periods, T\ and T2, depend sensitively on the dipole-dipole interaction and thus on the atomic 
distance r. They exhibit marked oscillations which decrease in amplitude when r increases. These oscillations seem to 
continue up to a distance of well over five wavelengths of the strong transition where we have stopped our simulation. 
The real part of the dipole-dipole coupling constant of the two V systems also exhibits damped oscillations. As a 
remarkable fact we find that these oscillations are in phase with those of T\ and T^. This correspondence is easy to 
understand intuitively since the real part of the dipole-dipole coupling constant enters the decay rates of the excited 
states of the combined two-atom system. 

In Section II we explain the methods employed which are based on the quantum jump approach [Q- |^| (equivalent 
to the Monte-Carlo wavefunction approach fl49| and to quantum trajectories ]50|). This approach is here adapted to 
two dipole-interacting V systems. 
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In Section III we define in more detail the three types of fluorescence periods of zero, single and double intensity. 
This involves an averaging procedure, both experimentally and theoretically. We then present the results of our 
simulations. 

In Section IV the numerical results are interpreted, and it is shown that one can associate three distinct subspaces 
of states of the coupled system with the three types of fluorescence periods. During each period the coupled system 
remains in the corresponding subspace. In Section V we discuss and summarize our results. 



II. THE QUANTUM JUMP APPROACH 



In the quantum jump approach |44|-ff8||, the time development of an atomic system is described by a conditional non- 
hermitian Hamiltonian, H con( i, which gives the time development between photon emissions, and by a reset operation 
which gives the state or density matrix right after an emission. For a general iV- level system these have been derived 
in Refs. The derivation is adapted here for a system consisting of two atoms. Slight modifications arise here 

since the field operator appears with different position arguments. 

We consider two atoms, each a V configuration as shown in Fig. 1, with the levels \2}i and |3)i (i = 1,2) 
and fixed at positions r,. We define operators Sfj (i — 1,2; j — 2,3) in the two-atom Hilbert space by Sy = \j)u(l\ 
and S^j = \l)u(j\. For simplicity we consider the case where the dipole moments of two atoms are the same, i.e. 
i(1|Xi|j')i = 2 (1 |j)2 = Dy. If Dij is real then the angle it forms with the line connecting the atoms is denoted 
by dj. In general §j is defined through cos 2 ■dj = |(Dij, r)| 2 /r 2 Dfj, where r = r2 ri. We assume the laser radiation 
normal to this line so that the lasers are in phase for both atoms. The two lasers are denoted by L2 and I/3. We 
take zero detuning and Ej,j(r, t) = Re [Ecy exp{— i(uijit — k • r)}] where hujji is the energy difference between level 1 
and level j. Making the usual rotating- wave approximation and going over to the interaction picture the interaction 
Hamiltonian becomes 

2 3 

i=l j=2 k,s 

+ H.C.]+J?L. (1) 



with the coupling constants 



and laser part 



/ \ 1/2 

9i ^ = ie \*ok?) ( D w> ek .*)' ( 2 ) 



^ = ?EI>{4+^}- (3) 

The Rabi frequencies of the lasers are Slj — (e/^Dy • Eoj for j = 2,3, and they are the same for both atoms. 
The operator Hi implicitly contains the dipole-dipole interaction of the two atoms, as seen from the conditional 
Hamiltonian H con d further below. In the Power-Zienau formulation, which we have used above, this interaction is due 
to photon exchange |lj] . 

Conditional Hamiltonian and waiting times. As explained in Refs. [^4|^8|, -ff C ond is obtained (in the interaction 
picture) from the short-time development under the condition of no emission, i.e. from the relation 

1 - J ffcondAi = (0 pft |f7 7 (A*,0)|0 ph ) (4) 

where the right-hand side is evaluated in second order perturbation theory for At intermediate between inverse optical 
frequencies and atomic decay times. In a similar way as for a single atom |^4|-^7j one obtains for the system of two 
three-level atoms |5l[] 

h- 3 



Hcond ~ 2i [^2 Aj ( S £> S V + S £j S ij, 



3=2 



+ Cj(S+S^ + S+S^) +H L (5) 
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with the r-dependent coupling constants 



' (1 - COS 2 



ikjir 
1 



(1 -3 cos 2 



(6) 



which contain the dipole-dipole interaction between the atoms. The dependence of Cj on r is maximal for -dj = ir/2 
(see Fig. 2). In the following we will assume for the Einstein coefficients and Rabi frequencies the relations 

Q 2 < ft 3 , n 2 < and A 2 r> . (7) 

Then A 2 and ReC2 can be neglected in H coa< n] we will also neglect ImC2 which is allowed if r is not small compared 
to A21, as seen from Fig. 2. 

Let \g), |e 2 ) and |e 3 ) denote the states where both atoms are in the ground state and the excited states |2) and |3), 
respectively, and let \sjk) be the symmetric and i\ajk) be the antisymmetric combinations of \j)\k) and |fc)|j). Then 
Eq. (||) becomes 

-ffcond = —. [A3 (|s 23 )(s 23 | + loasXaas]) + (^3 + C 3 )\s 13 )(s 13 \ + (A 3 - C 3 )\a 13 ) (a 13 \ + 2A 3 |e 3 )(e 3 | 
3 

+ {^2V2in 3 (\g)( Slj \ + |«i J -><e i |) +«n a (|«i3><«33| + |ai3)(a 23 |) +in 3 (|si 2 )<s 23 | - \a a ) (033 1) + H.c.} . (8) 

3=2 

Without lasers the conditional Hamiltonian is diagonal in this basis. 

Between emissions the atomic time development is given by U con d(t, 0) = exp {— iH con dt/h} which is non-unitary 
since -ff C ond is non-hermitian. The corresponding decrease in the norm of a vector is connected to the waiting time 
distribution for emission of a (next) photon. If at t = the initial atomic state is \ip) then the probability Pq to 
observe no photon by a broadband detector (over all space) is given by Q-^7|] 

P (t;m = \\U coad (t,0)\i>)\\ 2 , (9) 

and the probability density w\ of finding the first photon at time t is 

«>i(t;IV» = - j t Po(f,\i>)) ■ (10) 

For an initial density matrix instead of \ip) the expressions are analogous, with a trace instead of a norm squared in 
Eq. (|). 

According to Eqs. @ and (|l0|), A 3 ± ReC 3 describe the decay rates of \s\ 3 ) and |ai 3 ), respectively. From this 
the well-known fact follows that two atoms with dipole interaction can decay faster or slower than two independent 
atoms (superradiance and subradiance |2^]). ImC 3 corresponds to a level shift of \s\ 3 ) and — ImC 3 to a level shift of 
|di 3 ), caused by the interaction between the atoms. 

Reset matrix. Now we determine the reset operation which gives the state or density matrix right after a photon 
detection. Let the state of the combined system, atoms plus quantized radiation field, be given at time t by |0 p h) p (0 p h|, 
i.e. the atomic system is described by the density matrix p and there are no photons (recall that the laser field is 
treated classically). If at time t + At a photon is found by a non-absorptive measurement the combined system is in 
the state 

2P>CM* + A<, t)|0 ph ) P (0 ph \uj(t + At, <)F> (11) 

where 2P> = 1 — |0ph)lyi(0ph| is the projector onto the one or more photon space (since At is in the above range and 
thus small one could directly take the projector onto the one-photon space). The probability for this event is the 
trace over Eq. ([ll]) . For the state of the atomic system it is irrelevant whether the detected photon is absorbed or not 
(intuitively the photon travels away and no longer interacts with the atomic system). Hence after a photon detection 
at time t + At the non-normalized state of the atomic system alone, denoted by lZ(p)At, is given by the partial trace 
over the photon space, 

Tl{p)At = 

tr ph (jP>E//(t + At,t)\% h ) p{% h \U\{t + Ai,t)JP>) • (12) 
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We call n(p) the non- normalized reset state 46 1. Proceeding as in Refs. fh|[47]| an< ^ using perturbation theory one 
obtains fell 

K(p) = Re C 3 (S^pS+ + Sa 3 pS+ ) 

+A 3 (Si 3 pS+ + S; 3 pS+ 3 ) . (13) 

The normalized reset state is TZ(p) = n{p)/ttn{p). By Eq. ( [ll|) the normalization of n(p) is such that tryi72.(p)At is 
the probability for a photon detection at time t + At when the (normalized) state of the atomic system at time t is 
p. The laser field does not appear in the reset state, just as in the case of a single atom 46|47[], since its effect during 
the short time At is negligible. 

By a simple calculation one checks that Eq. ( |l3| ) can be written as 

n(p) = (A 3 + ReC* 3 ) R +P R f + + {A 3 - ReC* 3 ) R-pR f _ (14) 

where 

R + = (S u + 5 2 - 3 ) /V2 

= If)( s i3| + |si3><e 3 | + (|si 2 ) (s 2 3| - |ai2)(a23|)/V2 
R. = (Sr 3 - /V2 

= \g)(ais\ + |ai3>(e 3 j + (|si 2 )<a 23 | + |ai 2 )(s 23 |)/V2 . (15) 

If p is a pure state, p — \ip){i/i\ say, then R±pR\_ are also pure states. This decomposition of n{p) is advantageous 
for simulations of trajectories. As pointed out above, A 3 ± ReC 3 describe the decay rates of |si 3 ) and |ai 3 ) to \g). 
The state |e 3 ) can decay both to |si 3 ) and |ai 3 ), with respective decay rates A 3 ±ReC 3 . The decay rate of the states 
|s2 3 ) and |a 23 ) is A 3 and is the same as in the case of two independent atoms. 

Simulation of a single trajectory. Starting at t = with a pure state, the state develops according to U con d until 
the first emission at some time t\ , determined from wi in Eq. ( |l0| ) . Then the state is reset according to Eq. ( |i"3| ) to 
a new density matrix (which has to be normalized), and so on. 

The decomposition of n(p) in Eq. ( |l4| ) allows one, however, to work solely with pure states which is numerically 
much more efficient. One can start with a pure state develop it with U con d until t\ to the (non-normalized) 
\tl>(ti)), reset to one of the pure states R±\if>(ti)) /\\ ■ | with relative probabilities given by the factors A 3 ± ReC 3 
appearing in Eq. (|l4|), and so on. The waiting time distributions are not changed by this procedure. 



III. FLUORESCENCE JUMPS FOR TWO ATOMS 



For a single atom in a V configuration the existence of dark periods is due to two widely different time scales in the 
times between two subsequent photon emissions (cf., e.g., Refs. HgH®. The smaller time scale ' T '' is of the 
order of A 3 , while the larger time scale, T", is the inverse of the smallest eigenvalue of H con d/i for a single atom. 
One can pick a time, T, say, with f C T « T", and if the time between two subsequent photon emissions is longer 
than T one may then define this as a dark period. The mean duration of such dark periods is essentially independent 
of T, if chosen as above, and is given by T". When the waiting time between two photons is less than T the atom is 
said to be in a light period. The average intensity in a not too short light period is that of a driven two-level system, 
i.e. in our case the 1 and 3 levels. 

If one has two independent, noninteracting, atoms the combined fluorescence is just the sum of the individual 
contributions. When both atoms are in a dark period one has a dark period of the combined system. If only one atom 
is in a dark period one observes a fluorescence period with intensity of that of a single two-level atom (single-intensity 
period), and if both atoms are radiating one observes a double-intensity period. However, due to fluctuations in the 
emission times the latter two periods are not sharply defined if the atoms are so close to each other that one cannot 
determine from which atom a particular photon came. To distinguish the periods one has therefore to use an average 
photon intensity, obtained by means of an averaging time AT. This AT has to be large enough so that the photon 
intensity doubles when both atoms are not in a dark period. If, on the other hand, AT is chosen too large one 
may overlook short fluorescence periods, and one will see more seemingly direct transitions between dark periods and 
double intensity periods (double jumps). The analytic treatment of fluorescence of two independent atoms is easily 
obtained from the single-atom case pa . 

For two dipole-interacting atoms which are sufficiently close to each other the photons cannot be attributed to a 
particular atom, either, and one has to consider the two interacting atoms as a jointly radiating system. To be able to 
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differentiate between different fluorescence phases one again has to average the photon numbers over a time interval 
AT to arrive at an intensity, and in the following we present the results of numerical simulations which have been 
obtained by the methods explained in Section II. 

Atomic distances of a few wavelengths. Fig. 3 shows the number of emitted photons per time A^ 1 , averaged over 
a time interval AT = I9O/A3. If the atomic distance is larger than a third of a wavelength of the fast transition, 
r > 3A31, i.e. &3ir > 2, one can clearly distinguish three types of fluorescence periods, dark periods (0), single-intensity 
periods (1), and double-intensity periods (2). For two noninteracting atoms these would correspond to radiation of 
no atom, one atom, and two atoms, respectively. However, in the case of interaction and small distance the system of 
two atoms radiates as a whole, and in general one cannot attribute periods 1 and 2 to radiation of individual atoms 
as in the noninteracting case. 

The transition between the periods occurs rapidly but not instantaneously. The duration of the periods is long 
compared to the atomic time scale. From a sufficiently long trajectory one can obtain the average lengths of the 
periods, denoted by To, T\ and T2. The results are shown in Fig. 4. The Rabi frequency of laser 2, which drives 
the weak atomic 1-2 transitions, has been chosen in such a way that for independent atoms one has To = 2OOO/A3. 
As seen in Fig. 4, To is essentially independent of the atomic distance and thus of the dipole-dipole interaction. In 
contrast to this the two light periods are strongly distance dependent. In Fig. 4 (a), T2 varies between I2OO/A3 and 
26OO/A3 and a similar behavior is also seen in Fig. 4 (b). The curve for T\ resembles that for T2, except for smaller 
relative variation. 

There is an interesting correspondence between the r dependence of T\ and T2 with that of Re C3 . As seen from a 
comparison of Figs. 4 and 5 the variations with r of all three quantities seem to be in phase. For I111C3 the variation 
is out of phase, as seen from Eq. (||). Since ReC3 influences the decay rates of the two-atom systems, this in-phase 
behavior suggests that the variation in the lengths of periods 1 and 2 are due to an r dependence of the decay rates. 

Small atomic distances. Our simulations have shown furthermore that for r < jX 3 i the intensity in period 2 
decreases and no longer reaches that of two simultaneously radiating independent atoms. For very small distance only 
periods and 1 remain. The reason for this will be discussed below. 

IV. INTERPRETATION OF RESULTS 

For a single atom (in a V configuration as in Fig. 1) one can associate light and dark periods with certain atomic 
states and density matrices. During a dark period the atomic state rapidly approaches the eigenstate of Hcond with 
smallest imaginary part of the eigenvalue |44|j35| , F^I ■ This eigenstate is very close to |2), up to terms of order f^^b/f^ 
and ^2/^3 p2| . Thus, in a dark period, the atom can be regarded to be approximately in the state |2). During a 
light period the atom can be regarded to be in the equilibrium state (density matrix) of the 1-3 subsystem driven 
by O3, again up to terms of the above orders ]52|] . A jump from one fluorescence period to the other corresponds to 
a transition between these atomic states, and such a transition is caused by laser 2. With il 2 — and A2 = there 
would be no transitions. This correspondence clearly carries over to the three fluorescence periods of two independent 
atoms. 

We are now going to suggest a similar correspondence between fluorescence periods and states for two dipole- 
interacting atoms. Fig. 6 depicts the Dicke states of the two- atom system (see Eq. (||)). Dashed arrows indicate the 
weak driving by laser 2, solid arrows indicate strong driving by laser 3 and decay, respectively. Now, for Q 2 = 0, i.e. 
no dashed arrows, the states in Fig. 6 decompose into three non-connected subspaces, namely one spanned by |e2) 
and the two others spanned by the four inner and outer states, respectively: 

subspace : |e2) 

subspace 1 : \s 12 ), \ai 2 ), |s 23 ), |a 23 ) 
subspace 2 : \g), \s 13 ), \a 13 ), |e 3 ) 

If the two atoms are in state |e2) then each of them is in its dark state, and thus no photon can be emitted. For 
5^2 = 0, the time development in subspace 1 is exactly the same as that for a system of two noninteracting atoms in 
the same subspace of states. This can be seen directly from Eq. (||) and it is also physically obvious since two atoms 
can only interact via photon exchange if none of them is in the dark state |2). The photon rate for subspace 1 is 
therefore that of a single two-level atom with levels 1 and 3. Subspace 2 corresponds to the level scheme of a system 
of two dipole-dipole interacting two- level atoms (with levels 1 and 3) , as for example recently discussed in Ref . [^6| . 
The photon rate of this system is, in good approximation, twice that of a single two-level atom, provided r > j A31. 
For smaller atomic distance the photon rate rapidly decreases to zero, due to the increasing level shifts. 

If the two-atom system has been in one of the subspaces 0, 1, or 2, it will quickly approach the corresponding 
equilibrium state for driving and distance r. If fl 2 =f= the additional weak driving will, from time to time and 
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in analogy to a single V system, cause transitions between the three subspaces 0, 1, and 2, and each transition will 
correspond to a jump in the fluorescence. 

Thus the fluorescence periods 0, 1, and 2 should correspond to (the equilibrium states of) the subspaces 0, 1, and 
2, respectively. This is verified by the numerical evaluation in Fig. 7. In the lower part of Fig. 7 a particular 
realization of an intensity trajectory with alternating periods of fluorescence is plotted. The three upper curves show 
the populations of the three subspaces corresponding to this realization, obtained by the conditional Hamiltonian 
and reset matrix of Section II. The agreement between fluorescence periods 0, 1, 2, and subspaces 0, 1, 2 is striking. 
During dark periods the two-atom system is in the subspace 0, and similarly for periods 1 and 2. 

This correspondence, however, depends to some extent on how large AT is chosen for the averaging of photon 
counts. If AT is chosen too large one can overlook some jumps between subspaces and some very short fluorescence 
periods. If AT is chosen too small there may be large intensity fluctuations, resulting in an incorrect determination 
of the different periods. 

If the atomic distance decreases below ^31 the level shifts of |si 3 ) and |ai3) by ImC3 increase rapidly. This renders 
the driving by laser 3 within the subspace 2 much less efficient. Hence for very small atomic distances the driving 
is essentially restricted to subspace 1. This explains the vanishing of double- intensity periods for very small atomic 
distances. 



V. DISCUSSION AND SUMMARY 



Cooperative effects of driven two three-level atoms have been studied, where each individual atom can exhibit light 
and dark periods. The atoms were considered to be a fixed distance r apart. If r is of the order of a few wavelengths of 
the fluorescent light the individual photons are no longer attributable to a particular atom and the two-atom system 
radiates as a whole, due to the dipole-dipole interaction. In addition to dark periods the system shows two types 
of light periods, one with fluorescence intensity as if only a single atom were radiating, and the other with double 
intensity. 

We have proposed to study the mean duration, T\ and T 2 , of the two types of light periods as a quantity sensitively 
depending on the dipole-dipole interaction and thus on the atomic distance. Experimentally and numerically these 
quantities are easily accessible. 

We have performed fluorescence simulations for atomic separations of up to 5 wavelengths and have found oscillations 
in Ti and T2 of up to 40% in amplitude. The amplitude decreases with the atomic separation but the oscillations 
seem to continue for separations larger than 5 wavelengths. By simulations we have shown that the r dependence 
of Ti and T2 is in phase with that of the real part of the dipole-dipole constant. This is eminently reasonable since 
Re C3 directly influence the decay rates of the excited states of the two-atom system. 

We have associated the three types of fluorescence periods with certain subspaces of states for the two atoms and 
to equilibrium states in these subspaces. The equilibrium states depend on the driving of the strong transition and 
on the distance. The weak driving then causes transitions between the subspaces. The transition rates depend not 
only on the weak driving, but also on the form of the respective equilibrium states and thus on the strong driving and 
on the atomic distance. In contrast to the mean durations of the light periods the mean duration of the dark periods 
is practically independent of r. This is intuitively quite clear since in the dark state there is essentially no photon 
exchange and thus no induced dipole-dipole interaction. 

To define fluorescence periods one has to average the number of photon emissions over a time interval AT of a 
some finite length. Hence very short fluorescence periods are washed out and not observed, and this can lead to 
apparent direct transitions between double-intensity periods and dark periods, or vice versa, so-called double jumps. 
Experimentally these have been seen in Ref. ^7|. The cooperative effects of up to 40% found by us in the duration 
of the single-intensity and the double intensity periods are noticeable. Therefore we expect that also the frequency 
for the appearance of double jumps is modified by the dipole-dipole interaction, but we cannot predict whether the 
changes are of the two orders of magnitude reported in Ref. ^7j. Besides, the system in Ref. |27|] differs from the one 
considered here, and it is not obvious if and how our results would carry over to that system. 
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FIGURE CAPTIONS 

FIG. 1. V system with metastable level 2 and Einstein coefficient A3 for level 3. Vt 2 and f2 3 are the Rabi frequencies 
of the two lasers driving the weak 1-2 transition and the strong 1-3 transition, respectively. 



FIG. 2. Dependence of Cj/Aj (j = 2,3) on r 



FIG. 3. Number of photons, I(t), per time A3" 1 emitted by two dipole-interacting atoms, averaged over AT = 190/A 3 , 
for f2 2 = 0,01 A 3 , f2 3 = 0,5 A 3 , and ^3 = n/2. (a) k 3 \r = 10, (b) fc 3 ir = 5, and (c) k^r = 2. The dashed curve 
indicates the fluorescence type (zero, single, and double intensity). 



FIG. 4. T , 7i and T 2 as a function of r for i? 3 = tt/2 and AT = 250/A 3 . (a) Q 3 = 0.3 A 3 and (b) ft 3 = 0.6 A 3 . VL 2 
has been chosen such that for two noninteracting atoms one has T = 2000/A 3 . 



FIG. 5. ReC 3 /A 3 as a function of r for i9 3 = 7r/2. 



FIG. 6. Dicke states. The dashed and solid double arrows denote weak and strong driving, respectively. Simple 
arrows denote decays. 



FIG. 7 Correspondence between fluorescence types and subspaces 0, 1, and 2 for = 0,01 A3, SI3 = 0, 5A 3 and 
i9 3 = 7r/2. The atomic distance is chosen as 5A 3 i/27r. 



8 




9 




10 




10000 20000 30000 40000 50000 

t (units of I/A3) (a) 




10000 20000 30000 40000 50000 

t (units of I/A3) (c) 



11 



3000 
2800 
2600 
2400 
2200 
2000 
1800 
1600 
1400 
1200 
1000 



T -o- 
Tj -+-- 

T 2 -o- 



Vfe-J 

+ T' 



10 



15 k 31 r 20 (a) 



25 



30 





5000 




4500 










4000 


i — i 








o 


3500 










'3 


3000 












2500 




2000 



T -E3- 
Tj -+-- 
To O- 



10 



15 feir 20 (b) 



25 



30 



12 




13 




14 





population 


of subspace 










n 











population 


of subspace 1 



























population of subspace 2 








1 







I 1 1 1 1 1 1 

0.3 - - 

0.2 - I I I I I 

0.1 - 

1 L- 1 ^ — 1 ^ 1 "-L- 1 1 1 

20000 25000 30000 35000 40000 45000 50000 

t (units of 1/A 3 ) 



15 



16 



